Size fluctuations of the initial source and the event-by-event transverse momentum 
fluctuations in relativistic heavy-ion collisions* 
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We show that the event-by-event fluctuations of the transverse size of the initial source, which 
follow directly from the Glauber treatment of the earliest stage of relativistic heavy-ion collisions, 
cause, after hydrodynamic evolution, fluctuations of the transverse flow velocity at hadronic freeze- 
out. This in turn leads to event-by-event fluctuations of the average transverse momentum, (pr)- 
Simulations with GLISSANDO for the Glauber phase, followed by a realistic hydrodynamic evolution 
and statistical hadronization carried with THERMINATOR, lead to agreement with the RHIC data. In 
particular, the magnitude of the effect, its centrality dependence, and the weak dependence on the 
incident energy are properly reproduced. Our results show that bulk of the observed event-by-event 
(pr) fluctuations may be explained by the fluctuations of the size of the initial source. 
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We propose a new mechanism for generating the 
transverse-momentum fluctuations in relativistic heavy- 
ion collisions, based on the fluctuations of the initial size 
of the formed system and its subsequent hydrodynamic 
evolution. It is well established that a successful de- 
scription of the physics of relativistic heavy-ion collisions 
is achieved with the hel p o f relativistic hydrodynamics 
P, i, S i, i, i, 0, S S M HO, E3, effective at proper 
times ranging typically from about a fraction of a fm/c 
to a few fm/c, where statistical hadronization takes over. 
Numerous observables can be reproduced that way, such 
as the momentum spectra, elliptic flow, or the HBT cor- 
relation radii [HI, measuring the system's space-time 
extension. The initial condition for hydrodynamics is 
usually obtained from the Glauber approach, leading to 
the wounded-nucleon picture [13| (a wounded nucleon is 
a nucleon, that has collided inelastically at least once) 
or its variants [3, [l5[ ■ When the initial condition is ob- 
tained via Glauber Monte Carlo simulations, its shape 
fluctuates, simply reflecting the randomness in positions 
of the nuclcons in the colliding nuclei. 

In this paper we show that the event-by-event fluctu- 
ations of the initial size are substantial, even when we 
consider the class of events with a strictly fixed num- 
ber of the wounded nucleons, N w . The fluctuations 
are then carried over by hydrodynamics to the fluctu- 
ations of the transverse flow velocity at the hadronic 
freeze-out, which in turn generate the event-by-event 
fluctuations of the average transverse momentum, px, 
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FIG. 1: (Color online) Two sample collisions, both with equal 
number of wounded nucleons, N w — 198, displayed as circles. 
The left and right cases have (r) = 2.95 fm and 2.83 fm, 
respectively. 



of the produced hadrons. The mechanism is very sim- 
ple: a more squeezed initial condition leads to faster 
expansion, larger flow, and, consequently, higher (j>t), 
while, on the conrary, a more stretched initial condi- 
tion leads to slower expansion, lower flow, and lower 
(pr)- The event-by-event pr fluctuations have been a 
subject of intense theoretical and experimental studies 

m m m m m m m m, m m, m m m, m, m, 

as they may reveal important details of the dynamics of 
the system. Effects of inhomogeneities for various ob- 
servables have been studied in Ref. [3l|. The event-by- 
event fluctuations of the initial shape have been studied 
in detail for its elliptic component, where they lead to en- 
hanced elliptic flow p, m, h n m m m m, m sn . 

Our present study is similar in spirit, but it focuses on 
the size fluctuations. 

Consider the average transverse size, defined in each 
event as (we use the wounded nucleon model for the sim- 
plicity of notation) 



N w 

E 



(i) 



where x% and yi are coordinates of a wounded nucleon 
in the transverse plane. Examples of the spatial distri- 
butions of wounded nucleons are shown if Fig. [IJ where 
collisions of two 197 Au nuclei are viewed along the beam. 
We note that the two cases displayed in Fig. [T] although 
having equal numbers of wounded nucleons, have in- 
deed a rather different shape and size. The original 
positions of nucleons in each nucleus are randomly gen- 
erated from an appropriate Woods-Saxon distribution, 
with an additional constraint that the centers of nucle- 
ons in each nucleus cannot be closer than the expulsion 
radius d = 0.4 fm, which simulates the short-range re- 
pulsion. Nucleons from the two nuclei get wounded or 
undergo a binary collision when their centers pass closer 
to each other than the distance d — yj ijnn /tt, where 
cjnn is the inelastic nucleon- nucleon cross section. For 
the highest SPS, RHIC, and LHC energies it is equal to 
32, 42, and 63 mb, respectively. 

We introduce the notation ((.)) to indicate averaging- 
over the events. In order to focus on the relative size of 
the effect we use the scaled standard deviation, defined 
for a fixed value of N w as c((r))/ ((f)). The results of our 
Monte Carlo simulations performed with GLISSANDD are 
shown in Fi g. El We present the standard wounded nu- 
cleon model [1 31 ] (the three lower overlapping curves) and 
the mixed model [3,[II[ (three upper curves), where the 
Glauber source distribution is formed from the fractions 
(1 — a)/2 of the wounded nucleons and a of the binary 
collisions. The mixing parameter a is assumed to be 
equal to 0.12, 0.145, and 0.2 for the highest SPS, RHIC, 
and LHC energies, respectively. From the fact that the 
three curves for the wounded nucleon model overlap, we 
conclude that the effect is completely insensitive to the 
value of (Tnn within the considered range. The observed 
moderate dependence on the energy for the mixed model 
originates entirely from the different values of the a pa- 
rameter. For the binary collisions the size fluctuations are 
stronger than for the wounded nucleons, hence a larger 
value of a yield larger fluctuations. We note that the 
scaled variance of (r) is about 2.5-3.5% for central colli- 
sions, and grows towards the peripheral collisions approx- 
imately as 1/ \JN W . Very similar results are obtained for 
other variants of Glauber models, in particular models 
with overlaid distribution of particles produced by each 
wounded nucleon [l5[ . We have also checked that using a 
Gaussian wounding profile <JNN(b) [12] for the NN colli- 
sions, rather than the sharp wounding distance criterion 
applied here, leads to very similar curves. Furthermore, 
the use of the nucleon distributions including realistically 
the central NN correlations, as given in Ref. [43|], leads 
to practically no difference. In other words, the behav- 
ior displayed in Fig. [3] is robust, basically reflecting the 
statistical nature of the Glauber approach. 

The next step, crucial in converting the size fluctu- 
ation into momentum fluctuations, is hydrodynamics. 
We use the perfect boost-invariant hydrodynamics de- 
scribed in detail in Ref. [llj, which leads to a success- 
ful simultaneous description of the soft RHIC observ- 
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FIG. 2: (Color online) Event-by-event scaled standard devi- 
ation of the size parameter (r), evaluated at fixed values of 
the number of wounded nucleons, N w . The three lower over- 
lapping curves are for the wounded nucleon model at the NN 
cross section corresponding to the highest SPS (32 mb), RHIC 
(42 mb), and LHC (63 mb) energies. The three upper curves 
are for the mixed model at the subsequent energies. The 
dependence on the energy for the mixed model originates en- 
tirely from different values of the a mixing parameter, see the 
text for details. 



TABLE I: Shape parameters a and b, as well as the central 
temperature Ti, for various centrality classes for the highest 
RHIC energy of ^s^F = 200 GeV. 



c[%] 


0- 5 


5-10 


10-20 


20-30 


30-40 


40-50 


50-60 


60-70 


a [fm] 


2.70 


2.54 


2.38 


2.00 


1.77 


1.58 


1.40 


1.22 


b [fm] 


2.93 


2.85 


2.74 


2.59 


2.45 


2.31 


2.16 


2.02 


Ti [MeV] 


500 


491 


476 


455 


429 


398 


354 


279 



ables at midrapidity, such as the px spectra, elliptic flow, 
and the HBT radii, including the azimuthally-sensitive 
(azHBT) femtoscopy 



44] 



The essential ingredients of 
this approach are the Gaussian transverse initial condi- 
tion, early start of the evolution (to = 0.25 fm/c), and 
the state-of-the art equation of state [ll|, incorporat- 
ing a smooth crossover and interpolating between the 
lattice QCD results at high temperatures and the res- 
onance gas at low temperatures. The initial transverse 
energy-density profile is assumed to have the simple form 
e(x, y) = eo exp[— x 2 /{2a 2 ) — y 2 / (2b 2 )}, where x and y de- 
note the transverse coordinates. The width parameters a 
and b depend on centrality. To obtain realistic values we 
run GLISSANDO [l5| simulations, which include the eccen- 
tricity fluctuations. Then we match a and b to reproduce 
the values (a; 2 ) and (y 2 ) from the GLISSANDO Monte Carlo 
profiles. Thus, by construction, the spatial rms radii of 
the initial condition and its eccentricity are the same as 
from the Glauber calculation. The values of the used 
width parameters are given in Table H] The central value 
of the initial temperature is adjusted in such a way that 
the multiplicity of produced particles is reproduced. 
The average transverse momentum is also reproduced 



3 



by our model, as can be seen from Fig. [3l where we com- 
pare the model predictions (solid line) to the STAR data 
pa ], extrapolated to the full pr coverage. The dashed 
line corresponds to the model calculation with the STAR 
range 0.2 GeV < p T < 2 GeV used in Ref. 

Next, we analyze the hydrodynamic evolution with 
fluctuating initial conditions. Rather than doing an 
event-by-event calculation, which is tedious, it suffices 
(see the following) to see how much the predictions 
change when the a and b parameters are scaled by the 
values read off from Fig.O For instance, for c = 20—30%, 
which corresponds to N w ~ 165, we note from Fig.[2]that 
the scaled standard deviation of the size fluctuations is 
4.4% (for the wounded nucleon model). Thus we rescale 
a and b up and down by 4.4% and run the simulations. In 
addition, we also adjust the value of the central temper- 
ature Tj in such a way, that the energy contained in the 
profile is preserved. This is natural, as the total energy 
deposited in the transverse plane should be (up to possi- 
ble additional fluctuations) the same for a given number 
of elementary collisions. Hence a squeezed system has 
a higher central temperature than a stretched system. 
Thus, in some sense, we include also the temperature 
fluctuations discussed in Ref. [5f3]. Additional event-by- 
event energy fluctuations (in the considered pseudorapid- 
ity window I77I < 1) could be added on top of the analyzed 
effect, which would act as another source of momentum 
fluctuations, not included in this work. 

Hydrodynamics is run till the local temperature drops 
to Tf — 145 MeV jTl}, where freezeout occurs. The 
freezeout hypersurfaces in the space of the transverse ra- 
dius, r, and time, t, at z = 0, as well as the transverse 
velocities, for the 5% squeezed and stretched cases for 
c = 20 - 30% are shown in Fig. |U We note the fol- 
lowing features: the maximum expansion velocities, in- 
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FIG. 3: (Color online) Dependence of the average trans- 
verse momentum, {(pr)), on the number of wounded nucle- 
ons. The solid (dashed) line corresponds to the model calcu- 
lation in the full pr range (in the range of the STAR coverage 
0.2 GeV < pr < 2 GeV). The experimental points (obtained 
by extrapolating the pr distributions to the full range) are 
taken from Ref. [45(. 
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FIG. 4: (Color online) In-plane (solid) and out-of-plane 
(dashed) freezeout hypersurfaces in the transverse coordinate- 
time plane for c = 20 — 30%, with the indicated values of the 
maximum transverse flow velocity in units of c. Top (bottom): 
the 5% squeezed (stretched) initial condition. 



dicated by dots with labels, are about 10% larger for the 
squeezed case compared to the stretched case, in addition 
the squeezed case is somewhat more compact, i.e. the ra- 
dius r and the time t at freezeout are slightly smaller. 

The final stage of our calculation, turning the 
transverse-velocity fluctuations into the transverse- 
momentum fluctuations, is the simulation of the statis- 
tical hadronization at freezeout with THERMINAT0R [47| . 
which includes all resonances and decay channels from 
SHARE 48] . According to the Cooper- Frye formalism, the 
particles (stable and unstable, which subsequently de- 
cay) are formed at the freezeout hypersurface according 
to appropriate statistical distributions. For the squeezed 
case they are more "pushed" with the larger flow ve- 
locity than for the stretched case, thus they acquire a 
higher average transverse momentum, ((px)). The ex- 
perimental cuts of the STAR detector are used, with 
0.2 GeV < pt < 2 GeV. 
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In order to compare to the data, we analyze the STAR 



correlation measure [28 
shown to be equal to 



dyn 



(ApiApj). It can be 



7 dyn = C 2 Ot)) - "Tr - 
J * nv 



E 

fc=i 



(2) 



where fc is numbers the events, cf\(pt) the variance of 
Pt in a given event, and Nk denotes the multiplicity of 
the event. By construction, the second term cancels the 
uncorrelated (purely statistical) fluctuations in the first 
term, leaving in <Td yn the dynamical correlations only. 
Our simulations do not include the statistical event-by- 
event fluctuations of (pt), which would follow from the 
random nature of the statistical hadronization involving 
a finite number of particles. The procedure described 
below allows to avoid the tedious event-by-event stud- 
ies in the extraction of erdyn when it originates form the 
fluctuations of the initial condition. 

The full statistical distribution / ((pr)) in a given cen- 
trality class is a folding of the statistical distribution of 
(Pt) at a fixed initial size, centered around a certain Pt, 
with the distribution of Pt centered around ((pt))- The 
value of px, corresponding to a fixed initial size, fluctu- 
ates because of the fluctuations of the initial size. In the 
central regions both are close to gaussian distributions, 
hence we have to a very good approximation 



/«Pt» 



exp 



{(pi 



d p T x 



PT) 



(3) 



exp 



(Pt-((ptW 

2 ^dyn 



Carrying out the pt integration yields the distribution of 
(pt) centered around ((pt)} with the width parameter 



satisfying a = a. 



2 

stat 



dyn' 



With the above factorized 



form we may obtain the second term directly from the 



distribution of the initial size parameter 
bution is also approximately gaussian, 



/((r)) - exp 



2^«r» 



Its distri- 



(4) 



Because of the deterministic nature of hydrodynamics, 
Pt is a (complicated) function of (r) . Nevertheless, in 
the vicinity of the central values we have from the Taylor 
expansion 



p T - ((pr» = 



dpi 



d(r) 



«r) - «r») 



(5) 



(r)=((r» 



Substituting (J5]) into ([4]) and comparing to ([3]) we obtain 
the key result 



erdyn = <r((r» 



dpT 



d(r) 



;r)=«r)> 



or for the scaled standard deviation 
ffdyn <r((r)) «r)) dp T 



((p T )) «r» (( PT )) d(r) 



(6) 



(7) 



This result bears similarity to the formula derived by 01- 
litrault where er dyn /((p T )) — <r((s))/ ((s)) , with s de- 
noting the entropy density. The derivative in Eq. ([7]) can 
be computed numerically without difficulty by running 
two simulations at each centrality. We do it by compar- 
ing the average momenta obtained for the squeezed and 
stretched cases, as described above. Then 



o-dyn = ((pt))- - ((pt))- 
(( PT )) (( Pt ))- + ((pt))- 



(8) 



where — and + indicate the squeezed and stretched cases. 

Our final result is shown in Fig. where we com- 
pare the theoretical points to the experimental data 
from the STAR collaboration [28|. We note a strikingly 
good agreement between our calculation and the exper- 
iment, in particular for the standard wounded nucleon 
model. The mixed model, which is more realistic than the 
wounded-nucleon model, overshoots the data by about 
20%, producing even more fluctuations than needed. 
This may suggest that the coefficient dpT / d(r)\i r \—i/ r \\ 
in Eq. © is somewhat too large. The value of this 
coefficient incorporates all the dynamics (the initial con- 
dition, hydrodynamics, statistical hadronization). Mod- 
ifying these ingredients, not to mention the inclusion of 
viscosity effects, etc., may modify the value. Neverthe- 
less, we note a proper magnitude of the effect and the 
correct dependence on centrality, with an approximate 
scaling ad yn ((PT)) / ((pt)) ~ 1 / V N w . Also, since the re- 
sults of Fig. [5] very weakly depend on cr NN , to the extent 
that the hydrodynamic "pushing" is similar at various en- 
ergies, our results should weakly depend on the incident 
energy, which is a desired experimental feature. We re- 
mark that the described mechanism works independently 
of the charge of particles. In order to describe the charge 
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FIG. 5: (Color online) Comparison of the theoretical predic- 
tions for the scaled dynamical fluctuations, o-d yn / {{pt)) (for 
y^ijvjv = 200 GeV) , to the experimental data from the STAR 
collaboration [2^|. The lower (upper) crosses indicate our re- 
sults for the wounded nucleon model (mixed model). The 
experimental data range from y^ijvjv = 20 GeV (triangles), 
through 130 GeV (squares), 62 GeV (diamonds), to 200 GeV 
(dots). 
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dependence of fluctuations (e.g. as observed in Ref. [25[) 
other mechanisms are necessary. 

In conclusion, one can straightforwardly reproduce 
bulk of the dynamical event-by-event transverse momen- 
tum fluctuations, as measured at RHIC, with the mech- 
anism based on fluctuations of the initial size, which are 
then carried over by hydrodynamics to the fluctuations 
of the transverse flow velocity, and consequently to the 
transverse momentum of the produced particles. The 
hydrodynamic pushing is crucial in this scheme. With a 
realistic hydrodynamics, which has been earlier used to 
uniformly describe the soft RHIC data, our analysis in- 



dicates that there may be little room for other, truly 
dynamical, sources of fluctuations, such as (mini)jets 
[27l . [50l | or the formation of clusters at freezeout [5l|, H3] . 
Certainly, there are yet other sources of the pt correla- 
tions in addition to the afore-mentioned ones, such as the 
global momentum conservation, resonance decays, corre- 
lations from the elementary NN collisions in the corona, 
however, these should be considered at the "background" 
mechanism of the size fluctuations described in this pa- 
per. 
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